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Abstract — A relativistic 3-fluid 3D hydrodynamic model has been devel- 
oped for describing heavy-ion collisions at incident energies between few 
and ~ 200 A-GeV. In addition to two baryon-rich fluids which simulate 
mutually decelerating counterflows of target and projectile nucleons, the 
new model incorporates evolution of a third, retarded baryon-free fluid 
created by this decelerated baryonic matter. Different equations of state, 
including those with the deconfinement phase transition, are treated. A 
reasonable agreement with experiment is illustrated by proton rapidity 
spectra, their dependence on collision centrality and beam energy. 

1 Introduction 

During past twenty years, hydrodynamics proved to be quite a reasonable 
tool for describing heavy-ion collisions at moderate energies (see for example 
references in ^J). From general point of view this application of hydrodynamics 
seemed to be more successful at higher energies when the number of produced 
particles gets larger. However, as it has been noted by D.I. Blokhintsev J] in 
the beginning of the multiple-production era, the uncertainty principle intro- 
duces certain restrictions to this extension. Though a general conclusion of 
Ref . j2] , that the very beginning of the expansion stage should be considered 
quantum-mechanically, is valid till now, the modern quark-gluon picture of 
interactions shifts the results of Blokhintsev's simplified estimate to far ultra- 
relativistic energies. On other hand, it became clear that non-equilibrium 
processes are very important at this stage and that only a part of the total 
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collision energy can be thermalized. In this respect the standard one-fluid de- 
scription assuming instantaneous local equilibrium seems to be quite limited. 

In the present paper we demonstrate first results of a relativistic 3-fiuid 
3-dimensional hydrodynamic code developed for describing highly relativistic 
nucleus-nucleus collisions, i.e. up to energies reached at the CERN SPS where 
a partial deconfinement of hadrons may occur. Below, basic features of the 
hydrodynamic model will be presented and a numerical solution of 3-fiuid 
hydrodynamics will first be compared with observables. 

2 3-Fluid Hydrodynamic Model 

A specific feature of the dynamic 3-fiuid description is a finite stopping 
power resulting in a counter-streaming regime of leading baryon-rich mat- 
ter. This counter-streaming behavior is supported by experimental rapidity 
distributions in nucleus-nucleus collisions and simulated by introducing the 
multi-fluid concept. The basic idea of a 3-fiuid approximation to heavy- ion 
collisions Ollj is that at each space-time point x = (t,x) the distribution 
function of baryon-rich matter, /b(a;,p), can be represented by a sum of two 
distinct contributions 



initially associated with constituent nucleons of the projectile (p) and target 
(t) nuclei. In addition, newly produced particles, populating the mid-rapidity 
region, are associated with a fireball (f) fluid described by the distribution 
function ff{x, p). Note that both the baryon-rich and fireball fluids may consist 
of any type of hadrons and/or partons (quarks and gluons), rather then only 
nucleons and pions. However, here and below we suppress the species label at 
the distribution functions for the sake of transparency of the equations. 

Using the standard procedure for deriving hydrodynamic equations from 
the coupled set of relativistic Boltzmann equations with the above-introduced 
distribution functions fa (a =p, t, f), we arrive at equations for the baryon 
charge conservation 



( for a =p and t) and the energy-momentum conservation of the fluids 




(1) 




(2) 




(3) 
(4) 
(5) 




f;{x) + f:{x)-fi;{x)-fi^,{x). 
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Here = UaU'^ is the baryon current defined in terms of proper baryon density 
Ua and fiydro dynamic 4-veIocity normalized 1. Eq. (0) implies 

that there is no baryon-charge exchange between p- and t-fluids, as well as 
that the baryon current of the fireball fluid is identically zero, = 0. The 
energy-momentum tensors T^'^ in Eqs. (jSl)^© are affected by friction forces 
Fp , F^, and F^. Friction forces between baryon- rich fluids, F^ and F^, 
partially transform the collective incident energy into thermal excitation of 
these fluids (|Fp — F^\) and give rise to particle production into the fireball 
fluid (Fp + F^), see Eq. (jS)). F^^ and F^ are associated with friction of the 
fireball fluid with the p- and t-fluids, respectively. Note that Eqs. 
satisfy the total energy-momentum conservation 

d^iT^" + Tr + Tr) = 0. (6) 

In terms of the proper energy density, e^, and pressure. Pa, the energy- 
momentum tensors of the baryon-rich fluids {a =p and t) take the conventional 
hydrodjTiamic form 

Tr = {Sa + Pa) < < - g^'Po.. (7) 

For the fireball, however, only the thermalized part of the energy-momentum 
tensor is described by this hydrodynamic form 

Its evolution is defined by the Euler equation with a retarded source term 

^^yMMV(^) = j (a; - a;' - Uf{x')t) [f^{x') + F^{x') 

- Fi;{x)-Fax), (9) 

where r is the formation time, and 

^'^"^ ^ - \F,{x') + F,{x')\ 

is a free-streaming 4- velocity of the produced fireball matter. In fact, this is 
the velocity at the moment of production of the fireball matter. According to 
Eq. 0, the energy and momentum of this matter appear as a source in the 
Euler equation only later, at the time UpT after production, and in different 
space point x' — \Jf{x')t, as compared to the production point x'. 
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The residual part of T/^'' (the free-streaming one) is defined as 

The equation for T^^'^^^" can be easily obtained by taking the difference between 
Eqs. (0) and 0- If all the fireball matter turns out to be formed before freeze- 
out (which is the case, in fact), then this equation is not needed. Thus, the 
3-fiuid model introduced here contains both the original 2-fluid model with 
pion radiation OEIEI and the (2+l)-fluid model [71IH] as limiting cases for 
r oo and r = 0, respectively. 

The nucleon-nucleon cross sections at high energies are strongly forward- 
backward peaked. Since the involved 4-momentum transfer is small, the Boltz- 
mann collision term can be essentially simplified and in this case the friction 
forces, Fp and , are estimated as 

K = P,Pt [i< - O Dp + {u; + <) De] , (12) 

a =p and t, p =t and t =p. Here, pa denotes the scalar densities of the p- and 
t-fluids, 

Dp/E = rriN Kei (^P/Eispt), (13) 

where ttt-at is the nucleon mass, Spt = m% (up + Wt ) is the mean invariant 
energy squared of two colliding nucleons from the p- and t-fluids, V^gi = 
[spt(spt — 4m^)]^/^/2m^ is their mean relative velocity, and (Jp/Eispt) are de- 
termined in terms of nucleon-nucleon cross sections integrated with certain 
weights (see OElEl for details). 

Eqs. (I21)^(I1|) and ©, supplemented by a certain equation of state (EoS) and 
expressions for friction forces F^, form a full set of equations of the relativistic 
3-fluid hydrodynamic model. To make this set closed, we still need to define 
the friction of the fireball fluid with the p- and t-fiuids, Ff^ and Ff^, in terms 
of hydrodynamic quantities and some cross sections. 

To estimate the scale of the friction force between the fireball and baryon- 
rich fluids, similar to that done above for baryon-rich fluids, we consider a 
simplified system, where all baryon-rich fluids consist only of nucleons, as the 
most abundant component of these fluids, and the fireball fluid contains only 
pions. At incident energies from 10 (AGS) to 200 A-GeV (SPS) the relative 
nucleon-pion energies prove to be in the resonance range dominated by the 
A-resonance. The resonance-dominated interaction implies that the essential 
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process is absorption of a fireball pion by a p- or t-fiuid nucleon with formation 
of an i?-resonance (most probably A). As a consequence, only the loss term 
contributes to the kinetic equation for the fireball fiuid. After momentum 
integrating this collision term weighted with the 4- momentum p'^, we get ^ 

rp{eci)Ou 

Ff^(x) = Dfa p„ , (14) 

where transport coefficients take the form 

D,^ = iy^-^^(sfa)/(m;vm^) = ^^or^'M- (15) 

Here, V^^ = [(sfa — m% — ra^Y — Am^jm^Y^"^ / {2mMmT,) denotes the mean 
invariant relative velocity between the fireball and the a-fiuids, Sfa = {mT,U{ + 
rriNUa)'^, and cr^t^^is) is the parameterization of experimental pion-nucleon 
cross-sections. Thus, we have expressed the friction in terms of the fireball- 
fiuid energy-momentum density T°^, the scalar density pa of the a fiuid, and a 
transport coefficient Dfa- Note that this friction is zero until the fireball pions 
are formed, since Tj''^'^''^'^ = during the formation time r. 

In fact, the above treatment is an estimate of the friction terms rather than 
their strict derivation. As it is seen from Eq. (fT^ for the excited matter of 
baryon-rich fiuids, a great number of hadrons and/or deconfined quarks and 
gluons may contribute into this friction. Furthermore, these quantities may 
be modified by in-medium effects. In this respect, Dp/E and Dfa should be 
understood as quantities that give a scale of this interaction. 

3 Simulations of Nucleus-Nucleus Collisions 

The relativistic 3D code for the above described 3-fiuid model was con- 
structed by means of modifying the existing 2-fiuid 3D code of Refs. [31310] • 
In actual calculations we used the mixed-phase EoS developed in PUIITT] . This 
phenomenological EoS takes into account a possible deconfinement phase tran- 
sition of nuclear matter. The underlying assumption of this EoS is that un- 
bound quarks and gluons may coexist with hadrons in the nuclear environment. 
In accordance with lattice QCD data, the statistical mixed-phase model de- 
scribes the first-order deconfinement phase transition for pure gluon matter 
and crossover for that with quarks (TUIITT]. For details concerning the used 
EoS's, please, refer to [T2j. 

In Fig. 1 global dynamics of heavy-ion collisions is illustrated by the energy- 
density evolution of the baryon-rich fiuids [sh = + ^t) in the reaction plane 



6 V.D. Toneev, Yu.B. Ivanov, E.G. Nikonov, W. Norenberg, V.N. Russkikh 

of the Pb+Pb collision at i?iab = 158 A-GeV. Different stages of interaction at 
relativistic energies are clearly seen in this example: Two Lorentz-contracted 
nuclei (note the different scales along the x- and z-axes in Fig.l) start to 
interpenetrate through each other, reach a maximal energy density by the time 
~ 1.1 fm/c and then expand predominantly in longitudinal direction forming 
ge-like" freeze-out system. At this and lower incident energies the 
baryon-rich dynamics is not really disturbed by the fireball fluid and hence 
the cases r = and 1 fm/c turned to be indistinguishable in terms of ^b- 

Time-evolution of in Fig.l is calculated for the mixed phase model. 
Topologically results for EoS of pure hadronic and that of two-phase models 
look very similar. Due to essential softening the equation of state near the 
deconfinement phase transition, in the last case the system evolves noticeably 
slower what may have observable consequences jTUlE]- 

The energy released in the fireball fluids is of an order of magnitude smaller 
than that stored in baryon-rich fluids and depends on the formation time. At 
realistic values of the formation times, r ~ 1 fm/c, the effect of the interaction 
is substantially reduced. It happens because the fireball fluid starts to interact 
only near the end of the interpenetration stage. As a result, by the end of the 
collision process it looses only 10% of its available energy at i?iab = 158 A-GeV 
and 30%, at -Eiab = 10.5 A-GeV. Certainly, this effect should be observable 
in mesonic quantities, in particular, in such fine observables as directed and 
elliptic flows. The global baryonic quantities stay practically unchanged at 
finite r p. 

To calculate observables, hydrodynamic calculations should be stopped at 
some freeze-out point. In our model it is assumed that a fluid element decouples 
from the hydrodynamic regime, when its energy density e and densities in the 
eight surrounding cells become smaller than a fixed value Efr- A value Sfr = 
0.15 GeV/fm^ was used for this local freeze-out density which corresponds to 
the actual energy density of the freeze-out fluid element of ~ 0.12 GeV / fm? . 
To proceed to observable free hadron gas, the shock-like freeze-out ^3] is 
assumed, conserving energy and baryon charge. 

Proton rapidity spectra calculated for Au + Au collisions are presented in 
Fig. 2 for E^yeam =6, 8 and 10.5 A ■ GeV . One should note that hydrodynamics 
does not make difference between bounded into fragments and free nucleons. 
In the results presented the contribution from light fragments (d, t, ^He and 
"^He) has been subtracted using a simple coalescence model ^HI- This proce- 
dure allows one to reproduce reasonably evolution of the spectra shape with 
changing the impact parameter h for all the energies considered. Some discrep- 
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Figure 1: Time evolution of the energy density, — Sp -\- Et, for the baryon- 
rich fluids in the reaction plane {xz plane) for the Pb+Pb collision (-Eiab = 
158 A-GeV) at impact parameter 6 = 2 fm. Shades of gray represent different 
levels of £b as indicated at the right side of each panel. Numbers at this palette 
show the £b values (in GeV/fm^) at which the shades change. Arrows indicate 
the hydrodynamic velocities of the fluids. 



8 V.D. Toneev, Yu.B. Ivanov, E.G. Nikonov, W. Norenberg, V.N. Russkikh 



Au+AuILq) Hadron gas boS 




y-Ycm y-ycm y-y. 

Au+AuILq) Iwo phase model boS 




-10 1 -10 1 -10 1 

y-ycm y-ycm y-ycm 

Au+AulEfl) Mixed phase model EoS 





Eo=6 AGeV 


b=2fm 


E„=8 AGeV 


b=2fm 


E„=10 5 ASeV 

b=2fm 


60 


( / 

1 // 


_^ 6fm Y 


( / 

\4 

+ / 


ifm \ ' 
6fm A; . 


i. 6fm V 


id 




••»•;:■< '\ 

.,Jfm,,/o° \\v\ 
■•J>('°°" \\\\ 




• •*(:■;■::.■' ^ 


/ '/ ■ - AW 


20 




, llfni/^,- \\\ 




^llfm/^ = ° ' \\ \ 


















-10 1 -10 1 -10 1 

y-Ycm Y-Vcm y-Vcm 



Figure 2: Proton rapidity spectra from Au + Au collisions at three bombarding 
energies and different impact parameters 6 (given in fm) Three panels corre- 
spond to different equation of state. Experimental points are from ^13j 



Three-Fluid Simulation of Relativistic Heavy-Ion Collisions 9 



iO F 



Mixed Phase EoS 



Hadron Gas EoS 




\ S 



A 



A. 



Ylab 



Figure 3: Rapidity spectra of protons from Pb + Pb collisions for the mixed 
phase and hadron gas EoS. Solid and dashed lines are calculated for impact 
parameters 2 and 4 fm, respectively. Experimental points are from |16j . 



ancy observed for peripheral collisions near the target and projectile rapidity 
are mainly due to a not-subtacted contribution of heavier fragments. As seen, 
the rapidity spectra are only slightly sensitive to the EoS used. The same con- 
clusions can be drawn from Fig.3 for higher CERN SPS energies. One should 
notice that a general agreement here for central collisions is much better than 
that reached in any other 2-fluid hydrodynamic model. This originates mainly 
from larger energy-momentum transfer in (jl2j) rather than from the account 
for the third fireball fluid. 

4 Conclusions 

In this paper we have developed a 3-fiuid model for simulating heavy-ion 
collisions in the range of incident energies between few and about 200 A-GeV. 
In addition to two baryon-rich fluids, which constitute the 2-fiuid model, a 
delayed evolution of the produced baryon-free (fireball) fluid is incorporated. 
This delay is governed by a formation time, during which the fireball fluid 
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neither thermalizes nor interacts with the baryon-rich fluids. After the for- 
mation, it thermahzes and comes into interaction with the baryon-rich fluids. 
This interaction is estimated from elementary pion-nucleon cross-sections. Im- 
plementation of different EoS, including those with the deconfinement phase 
transition, may open great opportunities for analysis of collective effects in rela- 
tivistic heavy-ion collisions. Unfortunately, in spite of reasonable reproduction 
of observable proton rapidity spectra in the wide range of bombarding ener- 
gies and centrality parameters we are unable to favor any of considered EoS. 
Analysis of more delicate characteristics is needed. This work is in progress 
now. 
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